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Abstract 

The thermodynamic effects of local structure fluctuations in glassformers are analyzed in 
terms of energy basins and inter-basin hopping. Depending on the time-scale of measurement, one 
observes short-time thermodynamic properties related to a narrow set of basins, or equilibrium 
properties that include structural relaxation. The inter-basin hopping is manifested by the 
fluctuations of the pressure, internal energy, and other thermodynamic short time characteristics. 
Formulas relating the inter-basin fluctuations of pressure and internal energy to the differences 
between the long- and short-time susceptibilities are found. Based on obtained relations, we 
discuss the relative sensitivity of the structure to temperature and pressure. 

PACS number(s): 64.70.Pf 

1. The physical picture of a supercooled liquid 

Extensive study , has revealed common qualitative features in the cooling and supercooling 
behavior of glassformers but also substantial differences in the details of this behavior (see 
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••••••• • and references therein). Upon cooling and supercooling, both the viscosity T] and 

structural relaxation time increase exponentially (strong glassformers) or faster (fragile 
glassformers). In the current literature, this behavior is explained, following the ideas of 
Goldstein 12 and Stillinger (see 13 and references therein), by a special topography of the potential 
energy landscape in the part of the 3N-dimensional configuration space of a glassformer occupied 
by the amorphous state. Here, N is the macroscopically large number of particles in the system. 
The energy landscape picture is widely used in the literature 14 , 15 , 16 , 17 , 18 , 19 , to explain many static 
and dynamic properties of supercooled liquids. It is assumed that the energy landscape of a liquid 
may be described in terms of energy basins separated by energy or entropy barriers, so that each 
configuration belongs to one and only one basin. A more detailed description of the energy 
landscape can be found in the cited literature. 

At very low temperatures, the lifetime b of a finite-size system in a basin may be large 
compared to the particle vibration period x. Then, the ergodicity is broken at times t<0 b . A general 

20 21 22 

discussion of systems with broken ergodicity may be found in the literature , , . Above the 
glass transition temperature, this situation of broken ergodicity may be realized only in small 
systems. Due to the global nature of the energy landscape picture, the assumption of a long basin 
lifetime b is always violated in a sufficiently large system . An inter-basin step assumes a particle 
rearrangement in the system. The number 3N of dimensions of the configuration space, as well as 
the number of independently rearranging clusters, are proportional to the particle number N. Then, 
the frequency of inter-basin hopping is a global characteristic of the entire system, proportional to 
N. 

A complementary local picture of liquids and glasses describes a small cluster randomly 
chosen in the system. Basic ideas of this picture were long ago discussed in the free volume theory 
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, , , (for a recent discussion of locality and more recent references see, for example, ref. , ). 
A particle of a cluster vibrates in a cage made by surrounding particles with a vibration period x; 
the amplitude of these vibrations is small compared to the distance between particles. The 
vibration motion conserves particle arrangement in the cluster. At times t~0, another component of 
thermal motion results in changes in cluster structure; this component is termed rearrangement 
motion. The time is then a local characteristic, the average lifetime of a stable arrangement in a 
small cluster. We consider the liquid in the temperature and pressure ranges where 0»x. 

The glassformer behavior suggests that rearrangement changes in the local structure are 
thermally activated 29 , 30 . A fundamental suggestion is that global structural changes may be 
described in terms of local elementary events, thermally activated local rearrangements. Here, 
local means that while particle arrangement in a small rearranging cluster is significantly changed, 
particle displacements in the rest of the system, outside of the cluster, are small. 

The link between the global and local pictures is established when a local rearrangement is 
identified as an elementary step in the inter-basin thermal motion. Rearrangements represent a 
form of thermal motion in the liquid; the fluctuating characteristic in this thermal motion is the 
local structure. Fluctuations of the local structure have their manifestations in both kinetic and 
thermodynamic properties of the system. Below, we consider relations between structural 
fluctuations and thermodynamic properties at different time scales. For equilibrium systems, 
relations between thermodynamic fluctuations and susceptibilities are well known (see, for example 
the textbook 31 ). Those relations were used to estimate the contribution of structural fluctuations to 
susceptibilities of glassformer in conditions when this contribution is dominant (see, for example, 
32 , 33 ). Below, we discuss structural and vibration contributions to susceptibilities, and give a 
straightforward derivation of these relations based on definition of basin thermodynamics. We then 
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discuss an interpretation of these relations in terms of short and long time measurements. In the 
last section, we apply the derived relations to estimate the sensitivity of the local structure to 
changes of external parameters. 



2. Statistical mechanics of a glassformer in terms of basins 

In equilibrium 34 , the probabilities to find the system in a small element of the configuration 
space are given by the Gibbs ensemble. By summing up the probabilities of all configurations 
belonging to a basin a, one finds the probability W(a) to find the system in this basin. Here, a is 
the basin label. One may interpret 33 the statistical ensemble as representing a large body made of 
identical N-particle systems. In this interpretation, different parts of the body occupy different 
basins, and W(a) is the probability to find a part in the basin a. For an ergodic system, W(a) also 
describes the fraction of time a part spends in this basin. General theorems of statistical mechanics 
justify the use of the ensemble W eq (a) to calculate equilibrium properties, including fluctuation 
characteristics. 

In the Gibbs ensemble, the probability dW(q) to find a system of classical particles in an 
element Dq of the configuration space is 



F-H(q) 



dW(q) = e ksT Dq, 
F = -k B TlnZ(T,V), 
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H(q) 




Here, q are the particles coordinates, H(q) the microscopic Hamiltonian of the system, and 
F(T,V) the Helmholtz free energy. It is assumed that physically equivalent configurations, 
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differing only by permutations of identical particles, are counted as one configuration; technically, 
this may be accounted for by dividing the result of independent integration over all particle 
coordinates by the number of corresponding permutations. Equilibrium (or long time- scale) 
pressure P, internal energy U, compressibility (3, heat capacity C, and other thermodynamic 
characteristics are related to the thermodynamic potential F(T,V) by standard thermodynamic 
formulas. 

As assumed in the literature (see cited references, and discussion in the next section), results 
of a short-time measurement are represented by the thermodynamics of a system confined in a 
basin. This basin thermodynamic is defined 16 by the basin free energy F a 

F a (T,V) = -k B TlnZ a . 2 

H(q) 

Z a =je k ° T Dq, 3 

qea 

where a labels the basin. The basin partition function Z a is a sum over all intra-basin 
configurations qe a, and F a is the basin free energy. The probability W(a) to find the system in this 
basin is 

F-F a 

W(a) = e ksT , 

F a F 4 

Z = £ e kBT =e kBT 

a 

In (4), the partition function Z(T,V) is rewritten in terms of basin contributions; the sum goes over 
the set of all physically different basins. The basin free energy F a plays the role of the effective 
Hamiltonian 35 for the basin label a as a fluctuating variable. 
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One finds the basin thermodynamic characteristics by using the thermodynamic formulas 
with the basin free energy F a as a thermodynamic potential. For example, the basin thermal 
expansion coefficient Op for a system that occupies, at the time of measurement, the basin a (as 
explained in the next Section, this actually means a narrow basin set including a), is 

3 2 F, 

av 2 

A basin thermodynamic characteristic X(a) as a function of time fluctuates due to inter-basin 
(rearrangement) motion. At large times the system becomes ergodic. The large time average short- 
time value <X(a)> can be calculated as ensemble average 

<X(a)>=XX(a)W(a). 6 

a 

Both the average short time values and the equilibrium values of observables are 
characteristics of the thermodynamic equilibrium. 

3. Inter-basin fluctuations 

At time scales t»0 of a long-time experiment the liquid approaches thermodynamic 
equilibrium; the structural equilibration mechanism is inter-basin (rearrangement) motion. At 
equilibrium, the rearrangement motion may be described as an unbiased random flight of the 
system in the basin set {a}, with basin label a(t) as a random variable. The relaxation kinetics, and 
the correlation of fluctuations in time are determined by the distributions of both energy minima 
and barrier energies, while the equilibrium probability to find the system in a basin is universally 
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defined by the Gibbs distribution. Configurations with energies close to the top of the barrier have 
a negligible statistical weight at equilibrium, and correspondingly a negligible contribution to basin 
probability and to thermodynamics; the energy barriers only determine the time- scale at which a 
representative basin ensemble may be sampled. The universal Gibbs distribution for basin 
probabilities leads to universal relations between fluctuation characteristics and susceptibilities 
derived below. 

Inter-basin fluctuations change the local structure; the thermodynamic manifestations of 
these changes are fluctuations of P(t), U(t), and other short-time characteristics. In the statistical 

IT 

mechanics of equilibrium systems, there is a general scheme relating the fluctuation 
characteristics (binary correlation functions) to the susceptibilities (response functions). Below, we 
derive similar formulas that relate the magnitudes of inter-basin fluctuations to the differences 
between the average short-time and the equilibrium values of thermodynamic observables. 

Consider a basin thermodynamic quantity X(a) . The equilibrium value of X(a) describes 
the average of a long time (equilibrium) measurement of this quantity for a given part of the 
system. We use the notation X for this equilibrium value of X(a), and reserve the brackets <X> for 
the basin ensemble average 

F— F 

< X >= X X(a) W(a) =X X(a) . 7 

a a 

For pressure P and internal energy U, the equilibrium values coincide with the ensemble-averages 
of short-time values. For the pressure P, and basin pressure P(a), defined as P=-(3F/3V) T , P(a)=- 
(3F a /3V) T , one finds 

d — 
P = -^7(- k B Tln S ekBT ) = Z p ( a ) W ( a ) = < p ( a )) • 8 
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For the heat capacity c v , thermal expansion coefficient Op, isothermal compressibility (3t and the 
bulk modulus K=1/(3t, and other thermodynamic susceptibilities, the average short-time values 
differ from corresponding equilibrium values. The definitions, and measurements, of these 
quantities assume small changes in temperature and/or pressure. At a long time-scale, a change in 
temperature or pressure results in structural relaxation in the system, while the short time 
thermodynamics excludes contributions to susceptibilities coming from these structural changes. 
For the bulk modulus K= - V(aP/aV) T , K(a)= - V(aP(a)/aV) T , one obtains from (8) 

K = - V lf^rl w(a)--^£[P(a)-pfw(a). 9 



K is the long-time, equilibrium characteristic, while the first term on the right is the average short- 
time bulk modulus. Formula (9) may be rewritten as 



([£P(a)] 2 } = - — ^— [K- (K(a))] = - ^ AK . 10 

Here, we introduced the notation 8X(a) for the deviation of a basin characteristics from its basin 
ensemble average, and AX for the difference between the equilibrium value X and the average 
short-time value <X(a )> of the same characteristic: 

£X(a) = X(a)-<X(a) >; 

AX = X-^X(a)W(a) = X-(X(a)} • H 

a 

Similarly to the derivation of (9), one obtains for the thermal coefficient of pressure 
(3P/aT) v 



(£P(a)£U(a)) = k B T 2 A 



12 
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the internal energy U is defined by 



U = F+TS = F-T 



U(a) = F+TS(a) = F-T 



= (U(a)>; 
v dT L 
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By using the thermal expansion coefficient 0Cp=(1/V)(3V/3T)p related to (3P/3T) V by Kap= 
(3P/3T) V , K(a)ap(a)=(3P(a)/3T)v, one rewrites (12) as 

(<?P(a)£ U(a)) = k B T 2 A(^ p ) . 14 

The structural contribution to the thermal expansion coefficient was recently calculated by 
Stillinger and Debenedetti (ref.[16]). For the heat capacity Cv , one obtains 



C v =T 



as 

dT 



9F 

ydTj 



= <C v ) + ^([^U(a)] : 
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This relation may be rewritten as 



([^U(a)] 2 ) = k B T 2 VAc v , 16 

where c v = Cv/V. 

Formulas (10),(12),(14),(16) give the fluctuations-to-susceptibilities relations for inter-basin 
fluctuations. Thermal equilibrium at the long time scale is a necessary condition for validity of 
these relations. This condition may appear challenging for fragile glassformers in the vicinity of 
their glass transition because of a wide range of relaxation times in these systems . The fraction of 
the fragile glassformer that has the relaxation times larger than the time of long-time experiments 
increases upon supercooling. When this fraction become significant, one arrives at a crossover 
from short to long time behavior; when the fraction of structurally frozen clusters becomes 
dominant, only short-time regime may be realized. 
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The formulas derived in this Section relate independently measurable quantities, and thus 
provide a way to directly test the suggested mechanism of inter-basin fluctuations. Relatively small 
systems may be studied by the computer simulation method. New opportunities to experimentally 
study pressure fluctuations are provided by optical measurements involving chromophores (spectral 
hole dynamics, single molecule spectroscopy) in small systems imbedded in a rigid matrix. 
Indirectly, one may test predictions that are based on the derived relations. A prediction of relative 
sensitivity of the structure to changes in thermodynamic parameters is discussed in the next 
Section. 

4. Time scale for short-time kinetics. 

For a system with a long lifetime b »x in a basin, ergodicity at times t<0 b is broken. A 
general discussion of short time properties for such a system may be found in the cited ref. 20-22. 
Average value of the internal energy U, pressure P, volume V, and other quantities over a time t, 
x«t«0b, may be approximated as equilibrium properties for a system when it is confined in the 
basin. In a large system, however, the inter-basin fluctuations of pressure P(t) and other 
thermodynamic characteristics have a high frequency proportional to the macroscopic particle 
number N. We present here arguments that the basin thermodynamics may be applied to short-time 
measurements at times t<0; the time is a local characteristic independent from the particle number 
N in the system. 

At a short time-scale t«0, rearrangements take place in a small fraction ~ t/0 of the 
material, while in the rest of the material the structure is conserved. For additive thermodynamic 
quantities (volume, pressure, internal energy), contributions of different parts of the system are 
proportional to parts volume. Then, the relative contribution of those part where rearrangements 
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took place is proportional to the small volume fraction (~t/0«l) of this parts, and may be 
neglected. With this accuracy, one may describe the results of short-time measurement by the 
theoretical model of a system confined to a basin. The system actually changes basins during the 
short-time measurement, but the local structure for the basins visited by the system during the time 
of measurement coincide in most of the volume. The contribution of this structurally conserved 
part of the system to short-time thermodynamics coincides, up to small terms of the order of t/0, 
with the contribution of the system confined to any one of the visited basins. With this accuracy, 
the derived relations are between experimentally observable fluctuation and thermodynamic 
characteristics. 

5. Applicability to supercooled liquids 

It is well known that good glassformer may remain in a metastable supercooled state for a 
very long time, behaving at this stability time similar to stable systems. One may then approximate 
the basin ensemble for these supercooled liquids by the Gibbs ensemble. This approximation 
neglects the slow evolution of the metastable ensemble towards stable equilibrium state; this 

36 

evolution is considered in the nucleation theory . A metastable liquid remains in an 
approximately steady state only until the thermal inter-basin motion brings the system to a basin 
corresponding to crystalline clusters with the size about or larger than the critical nucleus size. 
Those states and basins will be referred to as crystalline. To model the metastable state, one 
considers a truncated configuration space that includes statistically significant structural states of 
the glassformer at the time when the metastable state is still relatively stable, and excludes 
crystalline states. Note that the probability to find an equilibrium liquid in a crystalline state is 
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negligibly small, so that the statistical mechanics for this system in the truncated configuration 
space coincides with that in the full configuration space. 

The formulas derived in above Sections are based on a Gibbs statistics. In the 
approximation of Gibbs statistics in the truncated basin set, these formulas may be applied to 
supercooled liquids. This application, however, assumes that equilibration in the truncated 
configuration space takes place at time scale when the system is steady. Changes in the metastable 
ensemble due to the growth of subcritical crystalline nuclei in the liquid matrix are small only at 
time-scales short compared to the lifetime of the metastable state. By using the approximation of 
Gibbs equilibrium in the truncated configuration space, we limit the applicability of our theory to 
supercooled states characterized by lifetimes much larger than the structural relaxation and 
equilibration times, so that the metastable system can sample a representative ensemble of local 
structures on the truncated set of basins. This sets a limit to the concept of a steady metastable state 
with properties independent from the thermal history of the system. 

6. Controlling parameters in supercooling 

The aim of this Section is to compare changes in the local structure caused by temperature 
and volume (or pressure) changes. Recent observations 37 , 4 have shown little difference between 
the increases of viscosity T| and local relaxation time upon cooling the liquid at constant pressure 
or constant volume. The conclusion of these studies is that the line 0,r|=const in the 
thermodynamic P,T plane is much steeper than the isochore V=const. 

Generally, any function Y(T,V) of the thermodynamic state defines a direction in the (T,V) 
thermodynamic plane of constant Y(T,V): the condition Y(T,V)=const in differential form is 
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3Y 

dV = p Y dT,p Y =-^|^. 17 

The parameter p Y determines the relative sensitivity of Y to volume and temperature changes: a 
temperature change dT results in the same change dY as the volume change dV=p Y dT. The 
pressure change dP Y =KdV/V that results in the isothermal volume change dV, equals 

dP Y =-— dV, 18 
V 

K= -V(3P/3V) t is the equilibrium bulk modulus. For example, the choice Y=V gives the isochoric 
pressure change dPv 

dP Y = Ka p dT , 19 
OCp is the thermal expansion coefficient. This relation gives a thermodynamic "compatibilizer" to 

compare changes in pressure and temperature. 

Changes in thermodynamic and kinetic properties reflect both changes at constant structure 

(changes in the vibration amplitude, density changes at constant structure (the mechanism better 

known in crystals), and structural changes. The general cause of structural changes is that, 

according to the general principles of statistical mechanics, the statistical weight of basins with 

lower free energies F a increases upon cooling; upon isothermal compression, the probability 

increases to find the system in a basin with more dense packing. The structure of the system is 

statistically described by the basin probability distribution W(a). When the temperature and 

pressure change, the probability W(a) to find the system in a basin a, and thus the structure, 

changes. For the relative change of the probability W(a) one finds 



8W(a) 1 



W(a) W(a) 



k B T 



- L ^ L ^dT+[P-P(a)]dV 
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For one basin, the condition SW(a)=0 determines the volume change dV a that compensates for the 
action of a temperature change dT: 

dV.=- U - U(a) dT. 21 
T(P- P(a)) 

The compensating volume change dV a depends on the basin chosen. One cannot expect to 
simultaneously compensate the probability changes 8W(a) for all basins. In other words, there are 
structure changes for any change of temperature and volume or pressure. For a given characteristic 
that depends on structure, for example for the viscosity r\, one can compensate the changes due to 
temperature increase dT by changes due to appropriate increase of pressure dP. The ratio dP/dT 
may be estimated by comparing structure changes caused by dT and dP. The ensemble average of 
the relative change 8W/W defined in formula (20) vanishes:<8W/W>=0. One considers then the 
ensemble average of (8W/W) 2 . A non- negative quadratic form Q, proportional to this ensemble 
average is 



Q(T,V) = X 



^ U(a) dT + ^P(a)dV n2 



W(a) . 22 



T 

Instead of volume change dV, one uses the corresponding isothermal pressure change dP=-KdV/V, 
to bring the quadratic form Q to the form 

Q(T, V) = A TT (dT) 2 -2 A TP dTdP+ App(dP) 2 , 23 

with 
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A TT =^(6U(a) 2 ); 



A T/ ,=-^(6U(a)6P(a)); . 24 
Ap,=^-((6P(a)) 2 } 

With the help of (10), (12), and (16), one writes the coefficients in terms of thermodynamic 
characteristics: 

A TP =-k B TVA(& P ),. 25 
K 



A DD =kJV 



IAKI 



l PP B ' • 

The condition Q=const determines an ellipse in the dP-dT plane. As a measure of 
compensation of structure changes due to dT and dP is given by the ratio 

s= (Apt) 26 
A A 

JT\. •ji'j" JT\. pp 

A complete compensation corresponds to s=l. A rather rough estimation of the coefficients in Q 
using extrapolation of the data on atactic polystyrene 38 , 39 , 40 at P~l bar, T=T g =150C gives s~ 1/500. 
This small value shows that the correlation of fluctuations of the energy U and pressure P is small, 
and one can not compensate structure changes caused by temperature change with those caused by 
a change of the pressure. 

A measure of the relative sensitivity of the structure to temperature and pressure is given by 
the ratio 



\A PP T I AKI 



from the mentioned data on polystyrene one gets 
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°K 



The coefficient q is larger but close to the isochoric coefficient (3P/3T) v =Ka P ~10bar/K. The 
accuracy of this estimation is low because the experimental data are for the glass transition range in 
the thermodynamic plane; as already mentioned, application of the theory to fragile glassformers 
close to the glass transition temperature imposes challenging conditions for the equilibration time. 
More accurate experimental data are needed to make quantitative predictions. 

1 3 

The slope of the glass transition line (line of constant viscosity r|=10 p) of polystyrene 
gives a substantially larger coefficient p n =(3r|/3T)/(3r|/3P) -40 bar/Kelvin. Some part of this 
coefficient may be explained by the fact that viscosity and relaxation time increase upon cooling 
even in a system with energy barriers not changing. At low temperatures, the data on viscosity of 
the liquid are traditionally fitted by the formula T|=r|oexp(A/kBT), with r| -0.01 poise and the 
effective excitation energy A(T,P) a function of temperature and pressure. Similar formulas may be 
written for relaxation time and other kinetic characteristics. At the glass transition, by convention, 
T|=10 13 poise, hence A/k B T g ~ 35. The constant viscosity line in the (P,T) thermodynamic plane 
yields the relation 




29 



or, for the coefficient pr. 



drj A dA 

_ 3T = T~3T 
" 3/7 3A 

dP dP 
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Strong glasses ' are characterized by small changes of A(T,V) upon cooling, while in 
fragile glassformers A(T,V) increases significantly. One possible cause of this increase is the 
increase of barriers without structural changes, due to the increase in density. Another source of 
fragility is due to the dependence of the excitation energies on local packing; the structural changes 
upon cooling towards more ordered packing (usually also increasing the density) result in a sharp 
increase in A. Some details of both mechanisms are studied recently by Sastry 41 . We suggest (see 
discussion in 42 ) that structural changes are the main cause of fragility in liquids studied in the cited 
papers 18 ' 4 . Then, as a first approximation, one may neglect excitation energy changes at constant 
structure. As shown above, any change in thermodynamic parameters T and P changes the 
structure, and compensation is impossible. One can, however, always compensate changes in 
viscosity, or in the apparent excitation energy A; this last condition define the coefficient p A 

3A dA 

One speculates that the sensitivity of the excitation energy to temperature and pressure is close to 
that of the structure, and pA~q, where q is defined in (27). From this formula and the definition 
m(T)=3logr|/3(Tg/T), m(T g ) being the Angell fragility characteristic 6,7 , one can rewrite (30) as 



mMk Q X 

P " =PA; mMk B T-A 

B g 



32 



Here, M=lnl0~2.3. At temperatures close to the glass transition range, for fragile glassformers 
mM<A/k B T g , and thus p n ~pA. When considered as a function of temperature and pressure, m 
decreases upon heating or lowering the pressure; for fragile glassformers, the denominator in 
formula (32) may become small, and correspondingly the coefficient pr|> p A as observed. 
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